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LECTURE NOTES FOR 18.102, SPRING 2009 



Lecture 24. Thursday, May 7: Completeness of Hermite basis 

Here is what I claim was done last time. Starting from the ground state for the 
harmonic oscillator 

(24.1) H = + x 2 , Hu = u , u = e-* 2 / 2 
and using the creation and annihilation operators 

(24.2) A=4- + x, C=-4-+x, AC-CA = 2ld, H = CA + ld 

dx dx 

I examined the higher eigcnfunctions: 

(24.3) Uj = C j u = pj(x)u (c), p(x) = 2 j x j + l.o.ts, Huj = (2j + l)uj 

and showed that these are orthogonal, Uj 1 j ^ and so when normalized 
give an orthonormal system in L 2 (M) : 

(24.4) t 



J 2j/ 2 (j!)57T3 

Now, what I want to show today, and not much more, is that the e 3 form an 
orthonormal basis of L 2 (M), meaning they are complete as an orthonormal sequence. 
There are various proofs of this, but the only 'simple' ones I know involve the Fourier 
inversion formula and I want to use the completeness to prove the Fourier inversion 
formula, so that will not do. Instead I want to use a version of Mchler's formula. I 
also tried to motivate this a bit last time. 

Namely, I suggested that to show the completeness of the 's it is enough to find 
a compact self-adjoint operator with these as eigcnfunctions and no null space. It 
is the last part which is tricky. The first part is easy. Remembering that all the ej 
are real, we can find an operator with the ej ;s as eigenfunctions with corresponding 
eigenvalues Xj > (say) by just defining 



oo oo « 

(24.5) Au{x) = ej)ej{x) = ^ Vj(» / 

a — n a — n 



ej(y)u(y). 



3=0 3=0 



For this to be an operator we need Xj — > as j — > oo, although for convergence we 
just need the Xj to be bounded. So, the problem with this is to show that A has 
no null space which of course is just the completeness of the e'j since (assuming- 
all the Xj are positive) 

(24.6) Au = <=> u±ejV j. 

Nevertheless, this is essentially what we will do. The idea is to write A as an 
integral operator and then work with that. I will take the Xj — w 3 where we [0, 1). 
The point is that we can find an explicit formula for 

oo 

(24.7) A w u = J^w^ ej(x)ej(y) = A(w,x,y). 

3=0 



I struggled a bit with this in class but did pretty much manage to do it. 
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To find A(w,x,y) we use some other things I did last time. First, I defined the 
Fourier transform and showed its basic propertyL 

(24.8) T : L\R) — C£,(R), T(u) = u, 

fi(0= / e-^u(x), sup\u\ < \\u\\ L i. 
Then I computed the Fourier transform of uq, namely 

(24.9) (^■« )(0 = V^«o(0- 



Now, we can use this formula, of if you like the argument to prove it, to show 
that 

(24.10) v = e- x2/4 => v = V^e-« 2 . 
Changing the names of the variables this just says 

(24.11) e-* 2 = / e txs - s2 / 4 ds. 

2V 7r Jm 

Now, again as I discussed last time, the definition of the Uj's can be rewritten 

(24.12) Uj (x) = { -± +x y e -* 2 / 2 = e* 2 l\-^e-* 2 . 

Plugging this into (24.11) and carrying out the derivatives - which is legitimate 
since the integral is so strongly convergent - gives 

2 /2 



(24.13) Uj (x) = ^— [ (-isy 



'/Us. 



Now we can use this formula twice on the sum on the left in (24.7) and insert 
the normalizations in (24.4) to find that 

(24.14) Y,^)eAv) = Y,-^rr- l { % e^-y-^-^dsdt. 

j=0 j = Q JK 

The crucial thing here is that we can sum the series to get an exponential, this 
allows us to finally conclude: 

Lemma 19. The identity (24.7) holds with 

1 ( 1 — w , x2 1 + w 2 



(24.15) A( W ,x,y)=^ 7 ==e^[- w — j (x + y) -^—^ (x - y ) 
Proof. Summing the series in (24.14) we find that 

e x 2 /2+y 2 /2 r 1 

(24.16) A(w,x,y) = ^ — / cxp(--w;st + isx + ity - s 2 /4 - t 2 /A)dsdt. 

^ 1 JR2 2 

Now, we can use the same formula as before for the Fourier transform of uq to 
evaluate these integrals explicitly. I think the clever way, better than what I did in 
lecture, is to change variables by setting 



(24.17) s = (S + T) /V2, t = (S - T)/V2 dsdt = dSdT, 

- l -wst + isx + ity - s 2 /A - t 2 /A = iS°^J - j(l + w)S 2 iT^J - 1(1 - w)T 2 . 
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The formula for the Fourier transform of exp(— x 2 ) can be used, after a change of 
variable, to conclude that 

exp(ib — — -(l + w)S jab = — r cxp( — — -) 



f ,.^x — y I,, , m 2v / 7r , (x — y) 2 . 

/ exp (iT— =»-- (l-w)T 2 )dT= * - exp - ^ y > . 

Jr V2 4 V(l - «0 2(1 -w) 



(24.18) 

/ o™(iT X ~ ^ _ - M - „,\T 2 \dT = 

Inserting these formulae back into (24.16) gives 

(24.19) A( W , x, y) = i exp f- ( " + ^ - ~ v) \ + X -+ V - 

1 7 ^ y> ^rVT^ V 2(1 + w) 2(1 -w) 2 2 

which after a little adjustment gives (24.15). □ 
Now, this explicit representation of A w as an integral operator allows us to show 
Proposition 31. For all real-valued f G L 2 (R), 

oo 

(24.20) EK^I'HI/lli- 

j'=i 

Proof. By definition of 

oo 

(24.21) ^|K e ,)| 2 = lim(.f,^.f) 
so (24.20) reduces to 

(24.22) lim(/,^/) = ||/||| 2 . 

uijl 

To prove (24.22) we will make our work on the integral operators rather simpler 
by assuming first that / e C° (R) is continuous and vanishes outside some bounded 
interval, f(x) = in |x| > i?. Then we can write out the L 2 inner product as a 
doulbe integral, which is a genuine (iterated) Riemann integral: 

(24.23) (f,A w f)=J J A(w,x,y)f(x)f(y)dydx. 

Here I have used the fact that / and A are real- valued. 

Look at the formula for A in (24.15). The first thing to notice is the factor 
(1 — w 2 ) — \ which blows up as w — > 1. On the other hand, the argument of the 
exponential has two terms, the first tends to as w — ► 1 and the second blows up, 
at least when i-j/^O. Given the signs, we see that 

if e > 0, X = {(x, y); \x\ < R, \y\ <R,\x-y\> e} then 

(24-24) sup|A(tu,a;,2/)| -» as w -» 1. 

x 

So, the part of the integral in (24.23) over \x — y\ > e tends to zero asun 1. 

So, look at the other part, where \x — y\ < e. By the (uniform) continuity of /, 
given S > there exits e > such that 

(24.25) \x-y\ < e \f(x)-f(y)\ < S. 
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Now we can divide (24.23) up into three pieces:- 

(24.26) (f,A w f)= f A(w,x,y)f(x)f(y)dydx 
Jsn{\x-y\>e} 

+ / A(w,x,y)(f(x) - f(y))f(y)dydx 

JSn{\x-y\<e} 

J Sn{\x-y\<e} 

where 5= [-R,R] 2 . 

Look now at the third integral in (24.26) since it is the important one. We can 

change variable of integration from x to t = \J jz^{x — y) and then this becomes 

f A{w 1 y + tJ l —^ 1 y)f{y) 2 dydt 1 

Jsn{\x-y\<t} \ l + w 



(24.27) A[w , y+t J l —^,y) 

V 1 + w 

= V^Tw) cxp {-WT^) {2y + tV ~ v)2 ) cxp ("?) ' 

Here y is bounded; the first exponential factor tends to 1 so it is straightforward to 
show that for any e > the third term in (24.26) tends to 

(24.28) ||/||| 2 as w -» 1 since J e"* 2/4 = 2y/n. 

Noting that A > the same sort of argument shows that the second term is 
bounded by a constant multiple of S. So this proves (24.22) (first choose 5 then e) 
and hence (24.20) under the assumption that / is continuous and vanishes outside 
some interval [— R, R]. 

However, the general case follows by continuity since such continuous functions 
vanishing outside compact sets are dense in L 2 (R) and both sides of (24.20) are 
continuous in / € L 2 (R). □ 

Now, (24.22) certainly implies that the ej form an orthonormal basis, which is 
what we wanted to show - but hard work! I did it really to remind you of how we 
did the Fourier series computation of the same sort and to suggest that you might 
like to compare the two arguments. 



